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Abstract. - We study the behavior of self-propelled nano- and micro-rods in three dimensions, 
confined between two parallel walls, by simulations and scaling arguments. Our simulations include 
thermal fluctuations and hydrodynamic interactions, which are both relevant for the dynamical 
behavior at nano- to micrometer length scales. In order to investigate the importance hydrody- 
namic interactions, we also perform Brownian-dynamics-like simulations. In both cases, we find 
that self-propelled rods display a strong surface excess in confined geometries. An analogy with 
semi-flexible polymers is employed to derive scaling laws for the dependence on the wall distance, 
the rod length, and the propulsive force. The simulation data confirm the scaling predictions. 
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Introduction. — Both in soft matter and in biol- 
ogy, there are numerous examples of swimmers and self- 
propelled particles. With a typical size in the range of 
a few nano- to several micro-meters, both low-Reynolds- 
number hydrodynamics [1] and thermal fluctuations are 
essential to determine their dynamics. Well-known bio- 
logical examples are sperm cells which are propelled by 
a snakc-likc motion of their tail [2], bacteria like E. coli 
which move forward by a rotational motion of their spiral- 
shaped flagclla [3] , and listeria which arc propelled by local 
actin polymerization at their surface [4,5]. In soft matter 
systems, synthetic self-propelled particles have been de- 
signed to perform directed motion. Examples arc bimetal- 
lic nanorods which are driven by different chemical reac- 
tions at the two types of surfaces [6-8] , or connected chains 
of magnetic colloidal particles on which a snake-like mo- 
tion is imposed by an external magnetic field [9]. 

Both in soft matter and in biological systems, surfaces 
and walls are ubiquitous. For example, bacteria in wet 
soil, near surfaces or in microfluidic devices [10, 11], or 
sperm in the female reproductive tract find themselves in 
strongly confined geometries. Already in 1963, Rothschild 
found that sperm accumulate at surfaces [12]. Thus sur- 
faces strongly affect the dynamics of swimmers and self- 
propelled particles. Typically, these particles live in an 
aqueous environment. Therefore, hydrodynamics plays an 
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important role in determining their behavior. The long- 
range hydrodynamic interactions (at distances from the 
wall much larger than the particle size, so that the particle 
can be approximated by a force dipole) induce a parallel 
orientation and effective attraction to the wall [13,14]. At 
short distances from the wall, the details of the propulsion 
mechanism become relevant. For example, it has been 
shown for E. coli that corkscrew motion of the flagclla 
leads to hydrodynamic attraction [15]. 

We study here the dynamics of self-propelled rod-like 
particles confined between two planar walls. Such particles 
capture the elongated geometry of most of the swimmers 
mentioned above. In the vicinity of a wall, the rod-like 
geometry of the particles is important, since it favors par- 
allel orientation — both with and without hydrodynamic 
interactions. We consider rods which are small enough for 
thermal fluctuations to play an important role. Thermal 
fluctuations induce a persistent-random-walk behavior of 
the trajectories in the bulk, and an entropic repulsion near 
the wall. 

In order to study these effects, we employ a particle- 
based mesoscalc hydrodynamics technique, in which hy- 
drodynamic interactions (HI) can be switched on and off 
easily. In the absence of hydrodynamic interactions, the 
effect of the fluid on the self-propelled particle corresponds 
to a Stokes friction and thermal fluctuations, as described 
by Brownian dynamics (BD). Our main result is that self- 
propelled rods accumulate at surfaces, both with Brow- 
nian dynamics and full hydrodynamics. Note that this 
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Figure 1: A rod is modeled by three filaments, interconnected 
by harmonic springs with spring lengths 4 and l c = ■s/^h- 

result is in contrast to rods pulled parallel to a surface by 
an external force [16,17]. The aggregation of self-propelled 
rods is found to depend on the rod length and on the 
strength of the propulsive force. We employ an analogy 
with scmi-flcxiblc polymers to derive scaling laws for the 
residence times near and far from the wall, and to distin- 
guish different dynamical regimes. 

Model and Simulation Technique. We model 
the rod of length / as a crane-like structure (see Fig. 1) 
[18]. Three semi-flexible filaments, each consisting of N m 
monomers, are arranged in a triangular cross section. The 
distances between the filaments and between monomers 
within a filament are the same bond length lb, so that 
I = (N rn — l)lb- The rod length is kept nearly constant 
by harmonic potentials Ub = if(|r^ i+ i| — lb) 2 for the bond 
vectors r^+i between neighboring monomers, and simi- 
larly for next-nearest neighbors. Spring constants K I 2 are 
chosen to be much larger than the thermal energy fcgT, 
so that length and bending fluctuations are irrelevant. 

To describe hydrodynamic behavior, we employ multi- 
particle collision dynamics (MPC), a particle-based meso- 
scopic simulation technique that naturally incorporates 
both hydrodynamic interactions and thermal fluctuations 
[19-22]. The MPC fluid consists of point- like fluid parti- 
cles of mass m in continuous space. The dynamics evolves 
in two steps. During the streaming step, particles prop- 
agate ballistically for a time interval h. In the collision 
step the particles are sorted into the cells of a cubic lattice 
with cell size a; they exchange momentum by a rotation 
of their velocities relative to the center-of-mass velocity 
v CTO of each cell by an angle a around a randomly cho- 
sen coordinate axis. We employ the parameters m = 1, 
k B T = l, a = l,h = 0.05, a = 130°, and pa 3 = 10 fluid 
particles per cell. This corresponds to measuring length 
and time in units of the cell size a and (ma 2 /hsT) 1 / 2 , 
respectively. Embedded particles can easily be coupled to 
the MPC fluid through inclusion of the monomers in the 
collision step. The crane-like structure of Fig. 1 with bond 
length lb = 0.5a and monomer mass M = 5m embedded 
in a MPC fluid represents a good approximation to slen- 
der rods in a Stokes fluid with an effective hydrodynamic 
radius of r = 0.45a [18]. With the fluid parameters given 
above, the viscosity is r\ ~ n^/rnksT/a 2 . This is large 
enough to keep the Reynolds number Re = plv/rj < 0.3 
for all considered rod lengths I and rod velocities v. 

An advantage of MPC is that hydrodynamic interac- 
tions can easily be switched off, while retaining similar 
thermal fluctuations and friction constants [23,24]. In this 



case, denoted random MPC, v cm is drawn from a Maxwell- 
Boltzmann distribution of variance kgT/p, see Ref. [24]. 
This allows us to separate Brownian dynamic from hydro- 
dynamic effects. 

Propulsion is achieved by applying a forward thrust 

Pi = for i > 1, (1) 

on each monomer i of the rod, where Pj^-i is the unit vec- 
tor connecting monomers i and i — 1, and ft is the force 
strength per monomer. When hydrodynamics is included, 
an appropriate reaction force is added to the fluid parti- 
cles in the collision cells occupied by monomers, such that 
the total momentum is conserved locally. This is equiva- 
lent to a rod that pushes the fluid at its surface towards its 
rear end. We employ a global thermostat to keep the tem- 
perature constant. For Brownian dynamics, the velocity 
v of the rod is proportional to /(. For our choice of pa- 
rameters, random MPC gives |u| ~ ft/lo with monomer 
friction coefficient 70 ~ 270 -^/mfc^T/a 2 . This leads to the 
Pcclet number 



with diffusion constant D = kgT / (Q^j l). For Pe 1, 
propulsion dominates and the rod should have a persis- 
tent directed motion, while for Pe ^ 1 thermal noise 
dominates and leads to a diffusive behavior. 

The rods motion is confined by two parallel hard walls 
a distance d apart. In the MPC simulations with hydro- 
dynamics, no-slip boundary conditions on the walls are 
implemented by a bounce-back rule (inversion of the ve- 
locity at the wall) and virtual wall particles [21]. In the 
directions parallel to the walls, we apply periodic bound- 
ary conditions. For random MPC, the bounce-back rule 
is applied to the monomers which hit the wall. 

Results. — The main result of our simulations is that 
self-propelled rods accumulate at a wall, both with and 
without hydrodynamic interactions, in contrast to passive 
rods which are depleted near a wall due to entropic rea- 
sons. We begin our study with the Brownian-dynamics- 
type random MPC. The rod accumulation at a wall is 
illustrated in Fig. 2, which shows the probability density 
P{z) to find the center of mass of a rod at distance z from 
the wall for various propulsive forces f t . Passive rods show 
a depletion layer of thickness 1/2; however, with increas- 
ing propulsion force /{, a pronounced peak develops near 
the wall. Note that this behavior emerges solely from self- 
propulsion - i.e. even without hydrodynamic interactions. 

Surface Excess. To quantify the surface localization, 
we define the surface excess 

fd/2 

8= / [P{z)-P b ]dz (3) 
Jo 

where Pb is the bulk probability density. A homogeneous 
distribution corresponds to s = 0, while full absorption at 
the wall corresponds to s = 1. For a passive rod, P{z) 
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Figure 2: (Color online) Probability density P(z) as function 
of the distance z from the surface, for various propelling forces 
ft. Simulations are performed with Brownian-dynamics-like 
random MPC. The corresponding surfaces excesses are s = 
— 0.11, s — 0.07, s = 0.21, and s = 0.33 for increasing f t . The 
rod length is I — 9.5a, the walls are located at z — 0a and 
z = 50a. A solid gray line shows the density profile of passive 
rods. 

increases linearly for < z < 1/2, and s — — l/(2d/Z — 1). 
To calculate the surface excess in the simulations, we de- 

1/2 

termine the probability p = J Q P(z)dz to find the center 
of mass of the rod within half the rod length from a wall. 
This can easily be converted into the surface excess via 

s = (dp-l)/(d-l) (4) 

with the assumption P& = {P{z))i/2<z<d/2i which is a 
good approximation since P(z) is nearly constant for 
1/2 < z < d/2, see Fig. 2. The surface excess is shown 
in Fig. 3 as a function of the rod length I for various 
propulsive forces ft- Passive rods with ft — show the 
expected negative surface excess due to entropic repulsion, 
in good agreement with the analytical equilibrium result. 
Long propelled rods show a strong surface excess, which 
decreases for smaller rod lengths (I /a < 10). Short and 
weakly-propelled rods behave like passive rods (with neg- 
ative s), but with increasing length, they show a crossover 
to a positive surface excess, with a minimum of the sur- 
face excess at intermediate rod lengths. This minimum 
occurs at a Peclet number Pe ~ 10. It should be noticed, 
however, that the dependence of the surface excess on the 
rod-length I and the propulsive force can not be combined 
into a dependence on the Peclet number alone — for ex- 
ample, because the curves for large f t in Fig. 3 depend on 
I but are essentially independent of ft- 

Scaling Arguments. In order to understand the mech- 
anism which is responsible for the effective surface ad- 
hesion of self-propelled rods, and to predict their behav- 
ior as a function of rod length, propulsive force, and wall 
separation, we exploit the analogy of the trajectories of 
self-propelled rods with the conformations of semi-flexible 
polymers. In the bulk, the rotational diffusion constant 




Figure 3: (Color online) Surface excess s as a function of scaled 
rod length I /a, for various propelling forces ft, as indicated. 
Simulations are performed with Brownian-dynamics-like ran- 
dom MPC. The (black) dashed-dotted line is the scaling result 
in the ballistic regime (see Eq. (8)), the (gray) dashed lines are 
scaling results in the diffusive regime (see Eq. (10)) for ft = 0.5 
and ft = 0.05. Wall distance is d — 50a. 

of a rod is D r ~ fc^T / (nl 3 ), which implies a persistence 
length 

C P ~ v/D r ~ nvl 3 /k B T (5) 

of the trajectory. The probability to find the self-propelled 
rod in a layer of thickness 1/2 near the wall can be ex- 
pressed as p — t w /{t w + Tb), where t w is the time the rod 
remains within this layer and Tb is the time it is located 
in the bulk (with 1/2 < z < d — 1/2). The surface excess 
is then obtained from Eq. (4). 

To estimate t w , we consider a rod, which at time t = 
is oriented parallel to the wall, and located very close to 
the wall with < z <C 1/2. As the rod moves forward, it is 
reflected when is hits the wall, and is thereby constrained 
to the positive half-space z > 0, see Fig. 4. This situation 
is very similar to a semi-flexible polymer, which is fixed 
at one end near the wall with tangent vector parallel to 
the wall; its bending rigidity k is determined by the persis- 
tence length, £ p = /t/feflT. In this case, the distance of the 
polymers from the wall increases as (z) ~ [ksT / n) 1 / 2 x i ^' 1 
and the orientation angle as (9) ~ (AbT/k) 1,,2 i 1 / 2 , where 
x = vt is the distance traveled parallel to the wall [25,26]. 
The condition (z) = 1/2 at t = t w then implies 




It is important to emphasize that there is no complete 
analogy between semi-flexible polymers and self-propelled 
rods. The main difference is that the equilibrium con- 
formations of long semi-flexible polymers touch a wall es- 
sentially tangentially - with a very small angle on both 
sides determined by the bending rigidity. The trajectory 
of a self-propelled rod, on the other hand, can approach a 
wall also almost perpendicularly, see Fig. 4. The rod will 
hit the wall and get stuck at first, but under the effect 
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the bulk for time Tb, either in the ballistic or in the diffusive 
regime. 



of a small wall slip and thermal fluctuations, it will then 
slowly reorient itself parallel to the wall. Thus, it should 
be noticed that we employ the polymer analogy only for 
those parts of the trajectory at the wall after the rod has 
oriented itself parallel to the wall. 

For the time Tb for the rod to stay in the bulk fluid, we 
have to distinguish two regimes. In the ballistic regime, 
with £ p 3> d, the rod travels essentially on a straight line 
between the walls, see Fig. 4. In this case, the bulk time is 



given by n 



d/ sin(6*), where 9 is the angle of the rod 



with the surface when it leaves the wall layer of thickness 
1/2. The polymer analogy explained above implies (9) ~ 
(7/£p) 1/3 for 9 < 1, so that 



d ftp 

V \ I 



1/3 



V 



knT 



1/3 



dl 2 ^v- 2 l\ 



(7) 



Thus, the scaling arguments predict in the ballistic regime 
the probability 

p = 1/(1 + a B d) (8) 

to find the rod in the wall layers, with a constant as which 
has to be determined numerically. Note that this expres- 
sion is independent of the velocity v, because both time 
scales t w and r& depend on v in the same way. 

In the diffusive regime, the persistence length £ p is small 
compared to the slit width d. In this case, the calculation 
of Tb is a mean first-passage time problem [27] . The mean 
first-passage time for a one-dimensional random walker 
starting at position zq to reach one of the walls is t± ~ 
zo(d— zq) To/^p, with the time unit To ~ £,p/v. For zq <C d, 
this implies Tb ~ zod. We take zq ~ £ p sin(0), because this 
is the distance from the wall layer which the rod travels 
before it forgets its orientation, and obtain 



Thus, in the diffusive regime, the scaling arguments imply 
P=- ..-,/.,. 7Z (10) 



1/3 



Tb 



I + a D d / t " 2/3 1- 4 / 3 



with a constant arj, where we have used that v ~ ft- Al- 
though p depends on ft (or v) in this case, this dependence 
decreases with increasing rod length I. 



Figure 5: (Color online) Rescaled wall time = r m ff and 
bulk time r 6 * = Tbff /10 for different propulsive forces ft, as 
indicated. Dashed lines are scaling predictions t w ~ l 5 ^ 3 and 
Tb ~ l 2 ^ 3 ■ The wall separation is d = 50a. 



In order to test our scaling predictions, we show in Fig. 5 

2/3 2/3 

the scaled times T w f t and Tbf t , as determined from the 
simulations. For sufficiently long rods, I > 5a, and not too 
small forces, f t > 0.5, the simulation results are found to 
agree very well with the scaling predictions (6), (7). We 
have also investigated numerically the d-dependence of Tb 
in the diffusive regime, and find good agreement with the 
linear dependence predicted by Eq. (9). 

The scaling results (8), (10) can now be used together 
with Eq. (4) to determine the surface excess s. With a sin- 
gle set of numerical prcfactors ajj = 0.23 and ao = 2.48, 
the Z-dependence in Fig. 3 can be described very well for 
various propulsion forces /(. This allows for the follow- 
ing interpretation. Simulations with strong propulsion, 
ft > 2.5, are in the ballistic regime. With weaker propul- 
sion, short rods are in the diffusive regime. With increas- 
ing rod length, the persistence length of the trajectory 
(Cp ~ v l 3 ) increases, leading to a crossover to the ballistic 
regime. For long rods, / > 20a, this leads to a surface 
excess independent of ft over more than two orders of 
magnitude, in excellent agreement with Eq. (8). 

Fig. 6 shows the dependence of the surface excess on 
the propulsive force f t in more detail. For small ft, the 
surface excess equals that of a passive rod. The surface 
excess increases, when the propulsive force becomes rele- 
vant compared to the thermal motion; it reaches a plateau 
when the persistence length becomes comparable to the 
wall separation d. These simulations results are again in 
good agreement with the scaling predictions. 

Hydrodynamic Interactions. In order to assess the ef- 
fect of hydrodynamic interactions, we compare in Figs. 6 
and 7 results of simulations with Brownian-type dynam- 
ics (BD) and full hydrodynamics. We focus here on rod 
lengths I < 10a in order to obtain good statistical accu- 
racy. The simulations show that the qualitative behav- 
ior is very similar as without hydrodynamic interactions. 
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Figure 6: (Color online) Surface excess s as a function of 
propulsive force f t . The parameters are I = 4.5 a, and d = 20 a 
or d = 50 a as indicated. The solid lines show the predicted 
asymptotic behaviors for small and large propulsive forces. The 
dashed lines indicate the scaling results in the diffusive regime 
(see Eq. (10)). The decrease of s for ft — 4 with HI is probably 
a finite- Reynold-number effect [Re ~ 0.3). 

However, hydrodynamic interactions reduce the surface ex- 
cess. This effect is most pronounced for intermediate rod 
lengths, see Fig. 7, and persists over a wide parameter 
range. 

For a rod moving far from a wall, hydrodynamic inter- 
actions increase the rotational diffusion coefficient and de- 
crease the translational friction coefficient logarithmically 
as a function of the rod length [28] , 

D r = K//2r) - 0.66 + 0{r/l)] , (11) 

7,1 = 27r^[ln(2/2r)- 0.21 + 0(r/Z)p 1 , (12) 

where 2r is the rod diameter. Thus a rod with hydrody- 
namic interactions (HI) should move faster than without 
at the same propulsion force, as it is indeed seen in our 
simulations. However, the persistence length of its trajec- 
tory, £ p ~ v/D r ~ f t /(j\\D r ), does not change, since the 
logarithmic factors cancel out (to leading order for long 
rods). When we insert these results in our scaling rela- 
tions, we find that both r b and reduced in the pres- 
ence of hydrodynamic interactions by a factor 1/ ln(Z/2r), 
while the ratio r w /r b remains unaffected in both the diffu- 
sive and the ballistic scaling regimes (to leading order for 
large I /a.) 

Fig. 8 shows our simulation results for the ratio 
t w (HI)/t w {BD) and T b (HI)/r b (BD) for systems with 
and without hydrodynamics. Clearly, hydrodynamic in- 
teractions speed up both processes. The observed decay 
with increasing rod length is in good qualitative agreement 
with the scaling prediction of a logarithmic dependence on 
the rod-length. 

However, the logarithmic factors in t w and r b cannot ex- 
plain the smaller surface excess with HI observed in Fig. 7, 
because s only depends on the ratio T w /r b . The reason for 



Figure 7: (Color online) Surface excess s as a function of scaled 
rod length l/a. Symbols show simulation results with (dashed 
line) and without (solid line) HI for different propulsive forces 
ft, as indicated. The dashed-dotted line shows the scaling re- 
sult in the ballistic regime (see Eq. (8)), the dashed line in the 
diffusive regime (see Eq. (10)). Wall distance is d = 20a. 

the decreased surface excess is that t w decreases more than 
r b , see Fig. 8. In order to better understand the behavior 
of these characteristic times, we show in Fig. 9 the distri- 
bution P of t w and Tf, for the same intermediate strength of 
the propulsive force as in Fig. 7. For Brownian-dynamics- 
type simulations, the distribution shows an initial decay 
at small r, followed by a second peak, both in t w and 
T/,. The two peaks in P{r b ) are caused by two different 
processes after the rod leaves the wall region, (i) reentry 
to the same wall after a short time under a flat angle, 
and (ii) crossing of the channel and reaching the other 
wall under a steeper angle. The same mechanisms are 
also responsible for the peaks in P(t w ), because a steeper 
entrance angle into the wall layer implies a deeper penetra- 
tion and thereby a longer residence time in the wall layer. 
The second peak gets larger when the rod gets longer and 
crossings are more likely than reentries due to the larger 
persistence lengths. 

For hydrodynamically interacting rods, Fig. 9 shows 
that the probability for short residence times in the wall 
and the bulk regions is significantly enhanced. No second 
peak is visible for rod length l/a = 7, but it develops for 
longer rods with l/a = 9.5. 

The reason for the stronger decrease of t w than r b when 
HI are switched on must be the hydrodynamic interac- 
tion with the wall. The far-field approximation of this 
interaction [13, 14] is certainly not quantitatively correct 
for the long rods and small wall separations we consider; 
however, it could still apply qualitatively. For 9 <C 1, 
it predicts a hydrodynamic torque proportional to p0/z 3 , 
with dipolc moment p ~ r/vl 2 , at a distance z from the 
wall [14]. This implies that as the rod moves towards the 
wall, it gets oriented parallel to the wall, and therefore 
penetrates less deeply into the wall layer. The torque 
can be expressed as an effective orientational potential 



J. Elgcti and G. Gomppcr 



1 i 1 1 1 1 1 1 r 




2 • 1 1 1 1 1 1 1 1 

23456789 10 

l/a 



Figure 8: (Color online) Ratios t w (HI)/t w (BD) and 
Tt(H I) / Tb(B D) of residence times for simulations with hydro- 
dynamic interaction and Brownian dynamics, as a function of 
the rod-length l/a. The propulsive force is ft = 0.5, the wall 
distance is d — 20a. The dashed-dotted line indicates the log- 
arithmic dependence expected from Eqs. (11) and (12). 

W r ~ p6 2 /z 3 . Thus, for small angles \0\ (and large dis- 
tances z) where W r < ksT, orientational diffusion still 
dominates. The combination of hydrodynamic alignment 
and small orientational fluctuations is consistent with the 
increase of short residence times in the distributions of 
Fig. 9 when hydrodynamic interactions are switched on. 

Summary. — We have shown that self-propelled rods 
in confined geometries show a strong surface excess. The 
surface excess is negative for a passive rod; it increases 
with the propulsive force, and saturates for large forces. 
The analogy with semi-flexible polymers allows the predic- 
tion of scaling laws, which are in good agreement with the 
simulation results. Our results for the aggregation of self- 
propelled rods at surfaces are relevant for many systems 
in biology and nanotcchnology. 
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